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Abstract 

We study an initial-boundary-value problem for time-dependent 
flows of heat-conducting viscous incompressible fluids in channel-like 
domains on a time interval (0, T). For the parabolic system with 
strong nonlinearities and including the artificial (the so called "do 
nothing") boundary conditions, we prove the local in time existence, 
global uniqueness and smoothness of the solution on a time interval 
(0,T*), where 0<T* <T. 



1 Introduction 
1.1 Preliminaries 

Let Q G C°'^ be a two-dimensional bounded domain with the boundary 
dQ. Let dQ = UT^ be such that F^i and F^r are open, not necessar- 
ily connected, the one- dimensional measure of F/j fl Fjv is zero and F/) 7^ 
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(riv = [J^rj^-*, r^Y* n r^^^ = for i 7^ j). in a physical sense, Q represents a 
"truncated" region of an unbounded channel system occupied by a moving 
heat-conducting viscous incompressible fluid, will denote the "lateral" 
surface and represents the open parts of the region Q. We assume that 
in/outfiow channel segments extend as straight pipes. All portions of Fat are 
taken to be flat and the boundary F^v and rigid boundary F^j form a right 
angle at each point where the boundary conditions change. Moreover, we 
assume that Td is smooth (of class C°°). 

The flow of a viscous incompressible heat-conducting constant-property 
fluid is governed by balance equations for linear momentum, mass and inter- 
nal energy [5] 

g {ut + {u ■ V)tt) - uAu + Vtt = ^(1 - a9)f, (1) 

div'u = 0, (2) 
CpQ {0t + u- V0) - kAO - ue{u) : e{u) = gaOf ■ u + h. (3) 

Here u = {ui, U2), tt and 9 denote the unknown velocity, pressure and temper- 
ature, respectively. Tensor e{u) denotes the symmetric part of the velocity 
gradient. Data of the problem follows: / is a body force and h a heat 

source term. Positive constant material coefficients represent the kinematic 
viscosity u, density g, heat conductivity k, specific heat at constant pressure 
Cp and thermal expansion coefficient of the fluid a. The energy balance equa- 
tion ([3]) takes into account the phenomena of the viscous energy dissipation 
and adiabatic heat effects. For rigorous derivation of the model like ([I])-© 
we refer the readers to pT] . 

Concerning the boundary conditions of the flow, it is a standard situation 
to prescribe the non-homogeneous Dirichlet boundary condition for tempera- 
ture and homogeneous no-slip boundary condition for velocity of the fluid on 
the flxed walls of the channel. In the case of temperature outflow boundary 
condition on Fat, we accept a frequently used assumption as "zero flux den- 
sity" , which is equivalent to the condition V6 ■ n = 0, sometimes referred to 
as a "do nothing" (or "natural" ) boundary condition and rather used in nu- 
merical simulations (cf. [21]). However, it is really not clear which boundary 
condition for u should be prescribed on Fat. The boundary condition 

du 

- TTTli + V- = FiUi (4) 

prescribed on F^-* is again often called "do nothing" (or "free outflow") 
boundary condition (cf. [7j, [8], [ID]). Here Fj are given functions of space 
and time and rij is the outer unit normal to ^ i = 1, . . . ,m. The bound- 
ary condition (jl]) results from a variational principle and does not have a real 



2 



physical meaning but is rather used in truncating large physical domains. It 
has been proven to be convenient in numerical modeling of parallel flows. 
For more information about application of this boundary condition and the 
physical meaning of the quantities Fi we refer to [7]. 

Remark 1.1 Assume that Fi are given smooth functions on , i = 1, . . . ,m, 

and consider the smooth extension F such that F{x,t)\^(i) = Fi(x,t). Intro- 
ducing the new variable V = it + F this amounts to solving the homogeneous 
"do nothing" boundary condition transferring the data from the right-hand 
side of (jl]) to the right-hand side of the linear momentum balance equation. 
Hence we will assume throughout this paper, without loss of generality, that 
= m (HD, i = 1, . . . ,m. 

Remark 1.2 To simplify the notation in the whole paper, we normalize ma- 
terial constants g, v, k, a and Cp to one. 

The paper is organized as follows. In Subsection 11.21 we introduce basic no- 
tations and some appropriate function spaces in order to precisely formulate 
our problem. In Section [2], we present the strong form of the model for the 
non-stationary motion of viscous incompressible heat-conducting fluids in a 
channel considered in our work, impose compatibility conditions on initial 
data, specify our smoothness assumptions on data and formulate the prob- 
lem in a variational setting. The main result of our work is established at 
the end of Section [2l In Section [3l we present basic results on the existence, 
uniqueness and energy estimates of the solution to an auxiliary linearized 
problems, the decoupled initial-boundary value problems for the Stokes and 
heat equations. The main result, stated in Section [21 is proved in Section H] 
via the Banach contraction principle. In the proof of local in time existence, 
presented in Subsection 14.11 we rely on the energy estimates for linear prob- 
lems, regularity of stationary solutions and interpolations-like inequalities. 
The global in time uniqueness of the strong solution is proved in Subsection 
14.21 using the technique of Gronwall lemma. 

Remark 1.3 Let us note that our results can be extended to the problem if 
we consider the so called "free surface" boundary condition onT^ and replace 
m by 

-Tin + v\Wu + {VuY]n = 0. 

However, to ensure the smoothness of the solution and exclude boundary sin- 
gularities near the points where the boundary conditions change their type 
some additional requirements on the geometry of the domain need to be in- 
troduced. This means that Vjq and Vd form an angle uj < n/A at each point 
where the boundary conditions change (see I23j ). 
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1.2 Basic notation and some function spaces 

Vector functions and operators acting on vector functions are denoted by bold- 
face letters. Unless specified otherwise, we use Einstein's summation conven- 
tion for indices running from 1 to 2. Throughout the paper, we will always 
use positive constants c, Ci, C2, . . . , which are not specified and which may 
differ from line to line. 

For an arbitrary r G [l,+oo], U{Q) denotes the usual Lebesgue space 
equipped with the norm || • ||L'-(n), and W'''P{fl), k > {k need not to be 
an integer, see [2D]), p G [1, +oo], denotes the usual Sobolev space with the 
norm || • ||vF'='P(n)- For simplicity we denote shortly W'^'P = W^'^{Qy and 



To simplify mathematical formulations we introduce the following nota- 
tions: 



In flSl)- ffTTl) all functions u, v, w, 6, (p are smooth enough, such that all inte- 
grals on the right-hand sides make sense. In ^ eij{u) denotes the compo- 
nents of the tensor e{u) defined by 





Let 



£u ■■= {u e C^inf; divn = 0, suppufirz) 



0} 



(12) 



and 



£e ■■= {0 e C°°{Q); supp^nT^^ 



0} 



(13) 
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and V^'P be the closure of in the norm of /c > and 1 < p < oo. 

Similarly, let Vq'^ be a closure of Eq in the norm of Then V^'^ and 

Vq'^ ^ respectively, are Banach spaces with the norms of the spaces 
and respectively. Note, that V^'^ 1//'^ V°'2 and respectively, 

are Hilbert spaces with scalar products (E]), O, (HU]) and (ITT]) , respectively. 
Further, define the spaces 

V^:={u\Se V°'2, au{u, v) = (/, v) for all ^ G Vi'^} (14) 

and 

:= I G V;°'', ae(^, ^) = (/i, ip)n for all G K/''} , (15) 
equipped with the norms 

and 

PWve ■■= \\h\\yO,., (17) 

where u and / are corresponding functions via dH]). Similarly, 6 and are 
corresponding functions via ( 1T5|) . 

The key embeddings W^'^ and Vg iy^'^(f2) are consequences 

of assumptions setting on the domain Q and the regularity results for the 
steady Stokes system in channel-like domains with "do-nothing" condition, 
see [21 Remark 2.2 and Corollary 2.3] and the "classical" regularity results 
for the stationary linear heat equation (the Poisson equation) with the mixed 
boundary conditions, see for instance [T9] . 

2 Formulation of the problem and the main 
result 

Let T G (0, oo) be fixed throughout the paper and Qt = f2 x (0, T), = 
Tjj X (0,T) and T^t = x (0,T). The strong formulation of our problem 
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is as follows: 



Uf -\- ylL • v J U — lAU, -\- V/ -r O J 


— J 


in n 




div u 


= U 


m Qt, 


(19) 


Ot + u-Ve - M - e{u) : e{u) 


= ef-u + h 


in Qt, 


(20) 


u 


= 


on Tdt, 


(21) 


e 


= 9 


on Tdt, 


(22) 


^ du 
-Vn + — 
on 


= 


on Ttvt, 


(23) 


06 
On 


= 


on Tnt, 


(24) 


uiO) 


= Uo 


in f2, 


(25) 


6(0) 


= eo 


in f2. 


(26) 



Here g is a given function representing the distribution of the temperature 6 
on Ft,, Uq and describe the initial velocity and temperature, respectively. 
Here we suppose that all functions in (fT5]) - fl26p are smooth enough and satisfy 
the compatibility conditions uq = and Oq = g on Td- 

Remark 2.1 Throughout the paper, fi denotes some fixed (arbitrarily small) 
positive real number (cf. fl27j) ). 

At this point we can formulate our problem. Suppose that 

/gL2+^(0,T;V°'2), heL'{0,T;V^''), (27) 
g e L\0,T;W'^'m, 9t e L'{0,T- L\n)), 
no e V.,, 9o E W^'^n), 9o - giO) E Ve. 

Find a pair [n, 9] such that 

niGL2(0,T;V°'2), n G ^^(0, T; I)„) n L°°(0, T; V^'^), (28) 
9t-gtE L2(0, T; V^^'^), 9 - g E L\Q, T; Ve) H 1^(0, T; V,''') (29) 

and the following system 

{ut,v) + au{u,v) + b{u,u,v) + {9f,v) = (/,«), (30) 
i9t,ip)n + ag{9,(f) + d{u,9,Lp) - e{u,u,Lp) - {9f ■ u,(f)Q = {h,(f)n{31) 

holds for every [v,ip] E V^'^ ^ ^^'^ and for almost every t E (0, T) and 

m(0) = no in Q, (32) 

^(0) = 9o in n. (33) 
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The pair [u,6] is called the strong solution to the system ([TS]) -f l2B]) . 

Let us briefly describe some difficulties we have to solve in our work. The 
equations (fT8|) - (l20|) represent the system with strong nonlinearities (quadratic 
growth of Vn in dissipative term e{u) : e(n)) without appropriate general 
existence and regularity theory. In [6], Frehse presented a simple example 
of discontinuous bounded weak solution U G fl H"^ of nonlinear elliptic 
system of the type AU = B{U, VU), where B is analytic and has quadratic 
growth in VU. However, for scalar problems, such existence and regularity 
theory is well developed (cf. [2T1 122]). 

Nevertheless, the main (open) problem of the system f|T8|) -f l26|) consists 
in the fact that, because of the boundary condition fl23l) . we cannot prove 
that b{u, u, u) = 0. Consequently, we are not able to show that the kinetic 
energy of the fluid is controlled by the data of the problem and the solutions 
of (|T8i) - (l26l) need not satisfy the energy inequality. This is due to the fact 
that some uncontrolled "backward flow" can take place at the open parts 
r^v of the domain Q and one is not able to prove global (in time) existence 
results. In [12]^ [H], Kracmar and Neustupa prescribed an additional condi- 
tion on the output (which bounds the kinetic energy of the backward flow) 
and formulated steady and evolutionary Navier-Stokes problems by means 
of appropriate variational inequalities. In [18] , Kucera and Skalak proved the 
local-in-time existence and uniqueness of a "weak" solution of the evolution 
Navier-Stokes problem for iso-thermal fluids, such that 

ut G L2(0,T*; Vi'2), uu G L\0,T*; (V^-^)*), < T* < T, (34) 

under some smoothness restrictions on no and V. In [25], the same authors 
established the similar results for the evolution Boussinesq approximations 
of the heat-conducting incompressible fluids. In [T7], Kucera supposed that 
the "do nothing" problem for the Navier-Stokes system is solvable in suitable 
function class with some given data (the initial velocity and the right hand 
side). The author proved that there exists a unique solution for data which 
are small perturbations of the original ones. 

In the present paper, we extend the results by Skalak and Kucera in 
|25j . However, our results are restricted to two dimensions. We shall prove 
local existence and global uniqueness of the strong solution to f|T8l) - fl26|) . i.e. 
for more general model than the Boussinesq approximations considered in 
cited references, and, moreover, such that i.a. [u,6] G L^(0, T*; W^'^) x 
L^(0,T*; ly^'^), which is strong in the sense that the solution possess second 
spatial derivatives. The main result of the paper is the following 
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Theorem 2.2 (Main result) Assume 

/gL2+'^(0,T;V°'2), heL\0,T;V,'''), (35) 
g E L'{0, T- W^^\n)), gt E L\0, T; L'{n)), 
UoEV^, 9oEW^'\n), 9o-g{0)EVg. 

Then there exists T* E (0,T] and the pair [u,6], 

ut E L\0,T;Y^/), u E L\0,T;V^) n L-(0, T; V^'^), (36) 
Ot -g^E L\0, T; V,'''), 6 - g E L^O, T; Ve) H L~(0, T; V^/'^), (37) 

such that [u, 6] is the strong solution of the system (ITHj) - (12^ . This solution 
is also globally unique. 

3 Auxiliary results 

Before we proceed to prove the main result of our paper, let us establish 
the well-possedness property for appropriate linear problems, which can be 
found in literature. 

Theorem 3.1 (The decoupled Stokes equations) Let f E L^{0,T;Yl'^) 
Then there exists the unique function u E L^{0,T]Vu) n L°°(0, T; V^'^), 
ut E L2(0,T; V°'2), such that 

{ut, v) + ttuiu, v) = {f , v) (38) 

holds for every v E V^'^ and for almost every t E (0,T) and u[0) = 0. 
Moreover 

II'"«IIl2{0,T;VS'') + ll**IU2(0,T;De) + II II L- {0,T;Vi'2) < c(f^) || / 1| i2(o,T;V2'') • (39) 

Theorem 3.2 (The decoupled heat equation) Let h E L'^{0,T;Vg''^) . 
Then there exists the unique function 9 E L\0,T;Ve) H L°°(0, T; V/'^), 
Ot E L2(0,T;1/g°'^), such that 

{et,ip)n + ae{e,ip) = {h,^)n (40) 

holds for every (f E Vg''^ and for almost every t E (0,T) and 6{0) = 0. 
Moreover 

Pt\\L^O,T;V°'^) + \\^\\L^O,T;Ve) + II ^ II L- (0,T;y;'2) ^ (^i^)\M L2(0,T;V°''')- (^1) 



8 



Remark to proofs of Theorem 13.21 and 13.11 The assertion of Theorem 
13.11 is proved in [31 Theorem 2.1] and [H Theorem 3.4]. We omit the proof 
of Theorem 13.21 since it can be established in the same way, see also [9l 
Chapter 5]. Note that the well-known approach to the proof is based on the 
Galerkin approximation with spectral basis and the uniform boundedness of 
approximate solutions in suitable spaces, see also [251 Chapter 3]. 



4 Proof of the main result 



4.1 Existence 

Let us introduce the following reflexive Banach spaces 



J\.rp 



{0|0eL2(O,T;P„), 0,gL2(O,T;VO'2), 0(0) = O} , (42) 
{^lj\^eL\0,T;Vg), ^teL\0,T-X''), V^(0) = O} , (43) 
{[0,^]|0GX^, V^ex^}, (44) 



respectively, with norms 



\m 



|[0,V']|| 



Xrj 



'tllL2(0,T;V2'^)' 



ll0l|L2{O,r;7?„) + 110^ 
||V'l|L2(0,T;1?e) + IIV^tllL2(0,T;y;'2), 



(45) 
(46) 
(47) 



Let us present some properties of X-^, Xf. and (consequently) Xt- Note that 
[21 Remark 2.2 and Corollary 2.3] yields the embedding T)u ^ W^'^, which 
implies 



X^ L~(0,r; W^'2) L°°(0,T;LP), 1 < j9 < oo. 
Let 4> ^ ^T- Using the interpolation inequality 

||0||w3/2.2 < C|| 

we get 



1/2 
Wl>2 



1/2 
W2.2 



(4^ 



(49) 



L4(0,T;W3/2,2) 



< C 



1/2 

L2(0,T;W2>2) 



1/2 

L°°(0,T;W1.2) 



< C 

where c = c(i7). Hence we have 

X^ L^(0, T; W^/'-') ^ L^(0, T; W^'^) L^(0, T; L^) 



(50) 
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for every 1 < p < oo. Similar properties hold for the space X^. Since [TU] 
yields the embedding Vq we get 

X| L~(0,T;iyi'2(l])) L°°(0,T;LP(fi)), 1 < p < oo, (51) 

and 

^ L\0, T- W^'^^\U)) ^ L^(0, T; lyi'^(n)) ^ L^(0, T; (52) 

for every 1 < p < oo. 

Let e L2(0,T; W^2,2(^^)) p L°°{0,T; W^'\n)). Raising and integrating 
the interpolation inequality 

\mt)\\w^+^-^in) < cii0(t)iii^2.2(f,)ii0(t)ii;;f;2(^), < J < 1, 

from to T we get 

mt)\\%..^n)'^ty < c ||0(t)||k^(^)||0(t)|lS^\:4Mty^ 

^ '^ll';^lll2(0,r;Vl/2.2(n))ll<^llL°°^(0,T;Vyl.2(Q)), (53) 

where c = c(r2). Hence we have 

L\0,T-W^'\Q)) n L^{0,T-W^'\Q)) ^ L^/^ {0,T-W^+^'\Q)) (54) 

^ L^/^{0,T-L°^{Q)) (55) 

for arbitrarily small positive j and consequently we deduce 

L'?(0,T;L°°) (56) 

and 

L'?(0,T;L°°(1])) (57) 

for 2 < g < oo. 

Remark 4.1 Setting [it>,i9] = [u — wq?^ — ^0^5'] ^^^s amounts to solving 
the problem with homogeneous boundary and initial data 

{wt, v) + au{w + Uq, v) + b{w + Uq,w + Uq, v) 

-if{'& + eo + g),v) = {f,v), 

(58) 

{-dt + gt, + a0{'d + 9o + g, ip) + d{w + Uo,d + 0q + g, ip) 

+ + 9o + g)f ■ {w + Uo),(p)n- e{w + Uo,w + uo,(p) = {h,(p)n, 

(59) 

w{0) = 0, (60) 
^9(0) = (61) 
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for every [v,(p] G V^'^ x Vg''^ and for almost every t G (0,T), where f G 
L2(0,T; h G L^{0,T; Vg^'^), G and 9o G P,. 

For arbitrary fixed [Wji}] G we now consider the nonlinear problem 

{wt,v) + au{w,v) = {f,v)-au{uo,v)-b{uo,Uo,v) 

—b{uQ, w, v) — b{w, Uq, v) — b{w, w, v) 
+(f(9o + g),v) + {fd,v), 

(62) 

{^t,v)n + ae{^,'-p) + e{w,w,ip) = -2e{w,Uo,ip) - e{uo,Uo,ip) + {h,ip)n 

-{gt, ip)n - aeiOo + g,(^)- d{uo, Oq + g, ip) 
-d{w, (f) - d{w, 00 + g,(p) - d{uo, ^, ^) 
+ ((^0 + g)f ■ uo, + (i^f ■ uo, 
+ ((^0 + g)f ■ w, (f)n + 0f ■ w, 

(63) 

w{0) = 0, (64) 
7?(0) = (65) 

for every [v, (f] G V^'^ x Vg''^ and for almost every t G (0, T). 

In the proof of the next lemma we verify that all terms on the right-hand 
sides of (162|) - (!63|) and the dissipative term e{w, w, ■) are well defined. From 
Theorem 13.11 and Theorem [32] we deduce that for arbitrary fixed [w,-t}] G Xt 
there exists [w,{}] = Z{[w,{}]) G Xt, the solution of the problem fl^ - flB^ . 
Consequently, the mapping Z : Xt Xt is well defined. 

Denote by Br{T) C Xt the closed ball 

Br{T) := {[0,^] G Xt; < R}- (66) 

We are going to show that the formal map Z : [w,i}] ^ [w,i)], where [w, i)] 
is a solution of the problem fl^^ - fl^S]l . has a fixed point in B^{T*) for T* 
sufficiently small. 

Let us prepare the existence proof. We now claim the following 

Lemma 4.2 For all R> 0, T > and for every [Wj-t}] G Bii{T) we have 

\\Z{[w,mx,<c^UT) + C2UT)uj{R), (67) 

where u is an increasing function depending solely on R, functions ^1 and C,2 
depend solely on T and ^i(T) — > 0+, C,2(T) — >■ 0+ forT 0+ and the positive 
constants Ci and C2 are independent ofT and R. 
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Proof The proof is rather technicaL Obviously, there exists a positive func- 
tion Ki(T) such that 

\\{f,-)-au{uQ,-) -b{Uo,Uo,-) + (/(^o + ^),-)|Il2(0,T;V2'2) 

+ Wih, ■)n - {gt, ■)n - ae{eo + g, ■) - d{uo,eo + g, ■) llL2(o,T;y;-2) 

+ ||e(no,no, ■) + iiOo + g)f ■ Uq, ■)n\\L2(o,T;V^''^) ^ cKi{T) (68) 

and KiiT) — )■ 0+ for T — )■ 0+, c is independent of T. 

Similarly and using the interpolation inequality we get 



||6(w;,Wo,-)IIl2(o,t;v2'^) < HIIl^II Vno||L4dt^ 



1/2 



< ||'"o||wi.4||'a7||L2(o_T;L4) 

< T^/''^||Uo||wi.4||'2^|U4(o,T;L4) 

< cT^'^u4j,Jw\\x^. (69) 
The other terms in fl62|l can be handled in the same way 

\\b{u, 

5 ""^5 ■)|Il2(o^T;V°'^) — ( / II'"'0|Il°°|| Viy||L2dt j 

< ||'"o||l°° ||iI'||l2(0,T;W1>2) 

< T^/^||'Uo||L-||'2^||L4(0,r;V\^i'2) 

< cT^'%u,yjw\\x^^, (70) 
||6(«;,'5;,-)|li2(o,r;v«-^) < II^IIl^II V5^||L4dt^ 

< r^''^||ii;||ioo(o_T;L4)||''i'||L4(o_r;Wi.4) 

< cT'"\\H\-, (71) 

and finally the last term in flB21) can be estimated using the interpolation 
inequality as 

ll(/^,-)IL2(o,T;Vr) ^ (^^ll/llL2|l^lli^(n)dt) ^ 

< r'^/^('+'^)('+'^)l||/||L2+.(0,T;L2)||^|| L2(4+n)/^(0,T;L°°(n)) 

< cT'^/[(2+'^)(^+^)]||/|U2..(o,r;L2)||^||x- 

(72) 
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The inequalities flB5]) - fl72]) and Theorem 13.11 yield the estimate 

\Mx^ < KiiT) + c(i?2 + R)T^'^ + cT^/I(2+M)(4+A<)]^^ (73) 

where Ki{T) — )■ 0+ for T — )■ 0+ and c is independent of T. Successively, we 
use (ISUj) and (175]) to estimate the dissipative term e(iu, lu, ■) in (IMl) 

l|e(w^,^^,-)|lL2(o,T;y;.2) < ||e(w)||^4dt^ 

^ II l|2 

^ II'"^IIl4(o,T;W1'4) 

< c||it>||^« 

< [Ki (T) + c(i?2 + i?)Ti/^ + cT'^/[(2+M)(4+M)] ^] 2 _ 

(74) 

Now we are going to estimate all remaining terms on the right-hand side of 
(163|) . We use interpolation inequality and (l50l) to estimate the dissipative 
term e{w, Uq, ip) 

l|e(w^,'"o,-)llL2(o,T;v;-2) < l|VWo|lL4||Vfi^||L4dt^ 

< T^/'^||VWo||L4||Vlt»||L4(o,r;L4) 

< cT^/%u4t,Jw\\x^ (75) 
and taking into account fl5UI) and flS2]) we arrive at the estimate 

IM(^,^,-)llL2(0,T;yr) < l|5^llL4||V?9||^4(f,)dt^ 

< T^/^ II lt> II ioo(o,T;L4) ll^ll L4(0,r;Vl/i.''(n)) 

< cTV4||iD|U^||^||^« 

< cTV^^(||^||^. + 11^11^.). (76) 
Let us proceed to estimate the remaining convective terms: 

IM(^,^0 + ^,-)|lL2(0,T;y;-2) < II^IIl4||^0 +^||^M(f^)dt^ 

< T^/^||ll>||ioo(o,T;L4)||6'o + fi'|U4(o,r;iyl.4(f7)) 

< cT'l^w\\x:^\% + g\\xf> 

< cTV^I (||{i;||^„ + ||^^o + ^7|li.) 

< cT'/'^Wwrx. + K.in (77) 
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where the positive function K2{T) — )■ 0+ for T — )• 0+. Further 
\\d{uo,^r)h^o,T-yr) ^ (/ \\^o\\l^\\Vnhd?j 

< cTVinoll^JI^II^.. (78) 
In the similar way one can deduce for the adiabatic terms the estimates 

< ||'*?||L4(4+M)/M(o,T;L°°{n))ll/llL2+M/2(o,T;L2)ll'"'llL4('*+A')/M(0,T;L°°) 
< C||??||^9 ||/||^2+M/2(o,T;L2)ll'»^l|x^ 

< cT^/t('+'^^^'-^"^^||^||x^ll/IU2+.(o,T;L2)||5;|U? (79) 

and, in the same way, 



% + g)f ■ W, ■)q\\^2(^0^T;V°''') 

^ 11^0 + 5'ilL4(4+M)/M(o,T;L°°(n))ll/llL2+M/2(o,T;L2)||'"^||L4(4+M)/M(o,T;L°°) 
< C\\0o + 5'llx|,ll/llL2+A'/2(0,T;L2(f7))ll''^l|x^ 

Finally, 



1/2 



||(^/-'"0,-)n|lL2(o,T;y;-2) < ll^lli-(n)ll/llL2|l'"0|lL-dt 

< cT^/'('+'^^^'+'^)^||^||x^ll/IU2+.(o,T;L2)||no||Loo. (81) 

Now the assertion is the straight consequence of inequalities fl68l) - fl8T|) . The- 
orem 13.21 and Theorem 13.11 The proof of Lemma 14.21 is complete. 

Lemma 4.3 ForT* G (0,T], sufficiently small, Z maps Bii{T*) into BuiT*) 
and realizes a contraction. 



Corollary 4.4 Lemma 4-2 , Lemma \4-3\ and the Banach fixed point the 



leorem 

yield the existence of a fixed-point [w,^] = Z{[w,'d]) in the ball Bji{T*) for 
sufficiently small T* e (0,7"]. By Remark \4-l\ the function [w, 0] = [w + 



Uq,^ + 9q + g] is the strong solution to the system f|T8l) - fl26|) on (0,T* 
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Proof of Lemma 14.31 Fix R > and choose r G (0,T] sufficiently small, 
such that (cf. dHZD) 

cMs) + C2Us)oo{R) <R Vs G (0, r]. 

Now we conclude that for every positive T* , < T* < t, Z maps Br{T*) 
into Br{T*). 

Now we prove contraction. Let ['2;i,'i9i], ["2^2; ''^2] ^ ^t, ^ = ^2 — Wi and 
a = '§2 — ^i- Theorem 13.21 and Theorem 13.11 imply the estimate 

\\Z{[w2,d2]) - Z{[Wi,di])\\xT <c(fi) (^116(2, Wo, ■)|lL2(o,T;Vr) 
+ IIK^O,2;,-)||i2(o,T;VS'^) + 11^(^2, 2;, •)|Il2(o,T;V«-2) + II ^(^> "^^l ' O II L2(0,T;VS'^) 
+ ll(/o^r)|lL2(o,T;Vr) + IM(^'^2,-)|lL2(o,T;y;-2) + . " ) || L2(o,T;yO'2) 

+ ||rf(z,^o ■)|Il2(o,T;V;-2) + IM(^0,C^, ■)|lL2(o,r;y;'2) 

+ ||e(w;2 -u?i,Wo,-)|lL2(o,T;y;-2)l|e(^2 -i'^i,^^2,-)|lL2(o,r;y;-2) 

+ ll((^0 + ^)/-2,-)n|lL2(o,T;y;-2) + ll(^/-'^l'-)f^llL2(0,T;y;'2) 

mf-z.-UW^io^T-y^))- (82) 
Following (jnED^dZHD we have (recall that c = c{Vi)) 



11^(2, -"0, 


■) L2(0,T;VS'^) 


< 


cT^/^llwolloJI^IIx^, 


(83) 


\\h{uQ,z, 


■) L2(o,T;V2'^) 


< 


cT^'^\\uq\\x,Az\\xj^, 


(84) 




■)l L2(0,T;VS'^) 


< 


cT^'^\\w2\\xj^\\z\\xi^, 


(85) 


116(2,-571, 


■) L2(0,T;VS'^) 


< 


cT^''^\\z\\xj^\\wi\\x;^, 


(86) 


II (/^, 


■)l L2(0,T;VJ1'2) 


< 


CT^/[(^^^^('^^)1||/IU2..(0,T;L2), 


lkllx«> (87) 


IM(2,^2, 


■) L2(o,T;yg"'^) 


< 


cT^/'Mx^^WUx^^, 


(88) 


\\d{Wi,(T, 


■)l L2(o,T;y°'2) 


< 


cT^'^WwiWxj^hWxf^, 


(89) 


11^(2,^0 


■) L2(o,T;yg"'2) 


< 


cT^I'\\z\\x^^\\e, + g\\xe^, 


(90) 


\\d{uQ,a, 


■)l L2(o,T;y°'2) 


< 


cT^'%Uq\\vA(7\\xO^. 


(91) 



Estimating the dissipative terms we arrive at 

||e(tU2 - it»i,'"o, ■)|lL2(o,r;y;'2) < c ||it»2 - lui ||x|; II^oUd^, (92) 
||e(i«;2 - ifi,it»2, •)llL2(o,T;yO'2) < c\\'W2 ~ Wi\\x;},\\w2\\x^, (93) 
||e(it»i,iU2 - -u^i, ■)|lL2(o,T;y;^2) < c\\'Wi\\xi^\\'W2 - Wi\\x^, . (94) 
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Following Theorem 13. II we can eliminate the term \\w2 — Wi\\x^ in fl^2|) - flM|l 
using the energy-like estimate (cf. dSH])) 

\\W2 - WiWx- < c{n) (^||6(z,no,-)|li2(o,T;vr) + IIK^^o, 2;, •)llL2(0,r;V2-') 

+ 11^(^2, ^,-)|Il2(o,T;V2'^) + ll^(2>l'-)llL2(0,T;Vr) + IK/^' O IIl2(o,T;VS'^)) 

(95) 

and then apply the inequalities (I83|) - (l87|) . Finally, the adiabatic terms can 
be estimated in the following way: 

II ((^0 + 5')/ ■ z, ■)n\\ i2(o,T;y;'2) 

< cT^/'^'+^^('+"^1||^o + ^7||x«II/IIl2+m(o,t;L2)||2||^., (97) 

||(f^/-'5^i,-)n||L2(o,T;y;'2) 

< cT'^/t^'+'^^^'+'^^^lk||^.||/|U2+.(o,T;L^)||5,i||x^ (98) 

and 

11(^^2/ ■ Z, ■)n|lL2(o,T;y;'2) 

< cT^/'('+^^^'+^^l||^2||x« II/I|l2+m(0,T;L2)||2||x- (99) 

Now let ['uji,'(9i], [w2,'&2] e Br{T) C Xt- Inequality (l82|) and the estimates 
(I83D-(I99D yield 

\\Z{[W2A]) - Z{[w^,d,])\\x, < cK{T)\\[w2,d2] - [w^A]\\xr. (100) 

where the positive function K(T) — )• 0+ for T — )• 0+ and c does not depend 
on T. Hence taking T* G (0,T] sufficiently small, the contraction of Z : 
Br{T*) — )■ Bii{T*) can be easily established. The proof of Lemma 14.31 is 
complete. 

4.2 Global uniqueness of the strong solutions 

Here we prove the global uniqueness of the strong solution stated in the 
main result. Suppose that all assumptions of Theorem 12.21 are satisfied and 
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there are two strong solutions [tii,6'i], [■U25 6'2] of f[TH]) -f l2B]) on (0,T). Denote 
z = Ui — U2 E and a = 9i — 62 E X^. Then z and a satisfy the equations 

{zt,v) + a^{z,v) + b{z,U2,v) + b{ui,z,v) - {fa,v) = 0,(101) 
{cTt, + + d{z, 9i, ip) + d{u2, cr, ip) - e{z, ui, v?) 

-e{u2,z,ip) - {af ■ ui,(p)n- {e2f ■ z,(p)n = 0(102) 

for every [v,Lp] G V^'^ x Vg''^ and for almost every t G (0,T) and z{0) = 
and cr(0) = 0. Hence substituting v = z and = cr we obtain estimates 

^^ll^WllviP + II^WIIvi.^ < |K^i(t),2(t),2;(t))| 

+ |Kz(t),n2(t),;^=(t))| + |(/(t)^(t),;2(t))| 

(103) 

and 



Id 
2dt 



n-JWimloa + IWrnl^.^ < \d{z{t),9,{t),am + \d{u2{t),a{t),am 



+ \e{z{t),ui{t),a{t))\ + \e{u2{t),z{t),am 

+ \{a{t)f{t)-u^{t),a{t))n\ 

+ \i92{t)f{t)-zit),a{t))n\ (104) 

for a.e. t G (0,T). To estimate term by term on the right-hand sides of 
(11031) and (I104p the Gagliardo-Nirenberg interpolation inequalities (cf. [U 
Theorem 5.8]) 

\\z(t)\W4 < c||zft)|L\/^||z^^^ii'/' 



L4 ^ i^\\zyi)\\^i,2\\zyi)\\j^2 , 

|k(^)l|L4(n) < C II (T(t) II J^i,2(f^) 

and the well-known Young's inequality with parameter 6 

ah < 5aF + C(5)6'? (a, 6 > 0, 5 > 0, 1 < g < 00, 1/p + l/q = 1) 

for C{5) = ((5p)^^/Pg~\ will be frequently used. 

Estimating the right-hand side of the inequality (I103p we arrive at 

\b{ui{t),z{t),z{t))\ < \\u,{t)U4z{t)\\^i,2\\z{t)U4 

< c\\u,mL4m\%U\m\\!" 
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\b{z{t),U2{t),Zm < \\u2{t)\M\z{t)\\l, 

< C\\u2{t) II W1.2 \\z{t) II W1.2 \\z{t) ||l2 

< 5||;2(t)||^..+C(5)||n2(t)||^,.||^(t)||L (106) 



and 



\if{t)ait),zm<\\fmi.4^mLHn)\\zit)h^ 

< s (iia(t)ii^..(^) + ii^(t)ii^..) + c{6)\\fmi (iia(t)iii.(^) + 



L2 



Now the estimates (fTOHj) and ffTnHjl - ffTOTD imply 
1 d 



(107) 



--Iiz(t)||^0,2 + ||2;(t)||^i,. < 6 (|k(t)||^,,.(f,) + ||2(t)||^i,2 

+ c{6) {wmwi + \\u,{t)\\i + \\u,{t)\\i,..,) \\zmi 

+ C(5)||/(t)||Llk(i)lli2(n). (108) 

Similarly, let us estimate all terms on the right-hand side of fll04p . Succes- 
sively 

\d{z{t),0,{t),crm < \\zmLAMimLAwmLHn) 

< c||^i(t) lk.2(^) ||a(t) 11^/^^^) ||a(t) ||^1,(^) ll^(t) 11^ 
<'5(lk(^)ll'vFi.2(^) + ||z(t)||^,,.) 

+ Cmemw..Hn){\WmlHn) + \\m\l) (109) 

and 

\d{u2{t),a{t),am < ||u2(t)||L4||Va(t)||L2||(T(t)|U4(^) 

< c||u2(t)||L4|k(t)||^L(^)||(T(t)||^(^(^) 

< 5|k(i)ll^.2(^)+C(5)||i.2(t)rL4|k(t)||i.(^). 

(110) 
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For the dissipative terms in (11041) we have 
\e{z{t),ui{t),a{t))\ < c||z(t)||wi.2|l'"i(^)llwi.4||o-(t)||L4(n) 

< c||z(t)||wi.2|l'"lWllwM||(T(t)||J/l2(^)||(T(t)||^i^(^) 

<5||z(t)||^.,+C'(5)||ni(t)||^.,.||a(t)||p^..(^)||a(t)|U.(n) 

(111) 

and in the same way we deduce the inequahty 

\e{u2{t),z{t),a{t))\ < c\\u2{t)\\w^A\\z{t)\\^i.2\\a{t)\\Li^a) 

< S (ll^WllwM + IkWII^M(n)) + C(5)1n2(t)||tvM|k(t)||i.(^). 

(112) 

The last two terms can be estimated as follows: 
Mt)fit)-u,it),ait))a\ < \\fit)h4MmL-hm%in) 

< c||/(t)||L2||ni(t)||L-||cx(t)||w/i.2(n)||a(t)||L2(f^) 

<'5|kWII^.2(f.) + cWII/WllLll^iWllLo.lk(t)lli.(f,) (113) 

and finally 

m{t)f{t) . zit),a{t))n\ < ||e2(t)||L^(n)||/(t)||LHI^WIlL.|k(t)IU4(n) 

<c\\e2mL^in)\\fmL4m\li,Um\]!^^^^^^ 



+ C(5)||e2(t)||io.(^)||/(t)||^.(||a(t)||i.(^) + ||z(t)||^.j. (114) 
Now the estimates (HM^ and (fT09|) - (fTTlD imply 
1 d „ 



+ lk(t)||J.^i,2 < S (llz(i)ll^i,. + ||a(t)f^,,.(^) 



+cis)[\m)\\ + \m)\\Un)\\fmi) 
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Choosing 6 sufficiently small and summing ( I108p and (IllSp we conclude that 

^ (ll^WIIvS'^ + MtWyO,) < X{t) {Wzml^, + Mt)\\lo,) , (116) 

where e L^{{0,T)). Now the uniqueness follows from the fact that 

z(0) = and a{0) = using the Gronwall lemma. 
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